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Abstract- The investment of connecting charging stations into 
distribution networks should be considered in the optimal 
planning of charging stations for electric vehicles (EVs). This 
paper proposes a mathematical model for multistage 
distribution network expansion planning taking optimal siting 
and sizing of EV charging stations into account. In the proposed 
model, the capacity of existing substations can be increased, 
existing network feeders can be replaced with wider section 
feeders, new feeders and new charging stations can be built. 
Network security constraints are considered. A new constraint is 
formulated to guarantee the topology of the obtained 
distribution network is radial. Numerical results indicate that 
the proposed method can allocate charging capacity among 
candidate charging station sites optimally, and the total 
investment and operating cost of the distribution network is 
reduced. 

Index Terms — Charging Station, Distribution Network, 
Electric Vehicle, Expansion Planning, Mixed integer linear 
programming. 



1. NOMENCLATURE 



A. Binary Variables 



x AK 
so 



Alternative RJ for replacing branch j at stage t 

Alternative AK for adding branch j at stage t 

Decision variable for adding or expanding 
substation / at stage t 
sl Capacity alternative SL for expanding substation / 

'•' at stage t 

^co Decision variable for adding charging station p at 

p,t stage t 



p '' p at stage t 

Vectors of decision variables at stage t and all 

xx 

' ' stages, respectively 

yf t Utilization of existing branch ; at stage t 

ro Utilization of alternative RO for replacing branch 



RJ 



j at stage t 

Utilization of alternative RJ for replacing branch 
j at stage t 

Utilization of alternative AK for adding branch j 
at stage t 

Vectors of branch utilization variables at stage t 
and all stages, respectively 

Flag whether node ; with load is connected at 
stage t 



B. Continuous Variables 

f F rF Current and maximum allowable current 

'•' '- /max ' i of existing branch ; at stage t 

rAK rAK Current of branch in the addition network 

' max * at stage and its maximum capacity AK 

Current of branch in the replacement 
f*>fmlj>fmLj network at stage t and maximum 
capacities of RO and RJ. 
Current vectors for the fixed, replaced, 
and added branches at stage t 
s s sl Power injection at node / at stage t and the 
Si,, , gmax,; > gmax,/ su b s tations ' maximum capacity. 

Column vector of the nodal injection at 
stage t 

Load shedding of node m and vector of 
load shedding at stage t, respectively 
Load of charging station p at stage t and 
the vector for all the charging stations 
Voltage magnitude of node /, m,p at stage 
t, respectively 
Column vectors of the nodal voltages at 
^f'^nin'^max stage t and their minimum and maximum 
limits 

c'" v (x t , P t , P t _ t ) Total investment at stage t 
c° per (r t ,y t ,P t ,) Total operational cost at stage t 

C. Parameters 



fF fR fK 
Jl >Jt >Jt 



8, 

r v 

P P 

P,l ' t 

V V V 

'/,»>' m,t> p,t 



expand Decision variable for expanding charging station p 



p.o 
p 

Demand ,t 



Whether node p has a charging station at 
initial stage (binary variable) 
The charging load at node p with a 
charging station 

Load of node m and all nodes at stage t 
Total charging demand at stage t 
Total number of planning stages 
Fixed expansion cost and variable 
investment cost for charging station p 
Present value factors for the investment 
S' t " v , S° per and operational costs at stage t 

II. Introduction 

There are growing concerns over energy security and 
greenhouse gas emissions. Transportation electrification is 
considered as an important way to solve these concerns. More 



' D 



and more electric vehicles (EVs) will therefore be integrated 
with electrical power systems through plug-in hybrid and 
electric vehicles. The integrations need charging facilities 
both at public and private places. The role of charging 
infrastructure is essential for the mass implementation of EVs 
as discussed in [1]. And the detailed requirements and 
conditions that should be considered for the design of 
charging infrastructure systems for EVs are also presented. 

Although charging at home is the first choice for many EV 
owners living in their houses, many other EV users who don't 
have private parking places or garages will have to use public 
charging facilities. Huge investment will be needed to build 
public charging places in countries like China [2]. A multi- 
objective planning model for the layout of electric vehicle 
charging station is proposed in [3], which takes into account 
factors including charging consumers' behaviour, distribution 
of the charging demand, factors of municipal planning. Case 
studies are carried out based on the background of Chengdu 
City in China. A method is given in [4] to optimally locate 
and size charging station for EVs based on grid partition. 
Genetic algorithm is employed to minimize the EV 
consumers' distance to the charging facilities. In [5], 
mathematical programming formulation is formed to 
determine the best locations for building alternative 
transportation fuel stations. EV charging is considered. 

To present, the work related to charging infrastructure 
planning mainly focuses on meeting charging demands with 
lowest investment costs and best users' convenience. The 
planning of charging locations doesn't consider the costs of 
connecting charging facilities to the grids. In this paper, the 
distribution network reinforcement planning is carried out 
simultaneously with the optimal siting and sizing of EV 
charging stations. 

III. Mathematical Formulation 

A. Problem Description 

The distribution network expansion planning problem 
considered in this paper takes the following factors into 
account [6, 7]: 

(1) the multistage planning problem has T stages with 
given duration; 

(2) the distribution network includes existing feeders 
connecting to the upstream substation(s), candidate branches, 
existing and new load nodes, candidate and even existing 
nodes for charging station(s); 

(3) to satisfy traditional load demands and future charging 
requirements, existing substation and charging station can be 
expanded, new substation and charging station can be added, 
and transfer capability can be increased by adding new lines 
or replacing conductors if necessary at any stage; 

(4) the operational and maintenance costs are calculated 
with each branch. 

B. Objective Function 

The objective function for the integrated distribution 
network and charging station expansion planning includes 



two parts: the investment cost c™ (x t , P t ,P t x ) and the 
operational cost c° per (r t , y t , P t , ) . The objective function is 

as follows: 

min C{x,y,r,P) = 

T (I) 

^rc'rix^PM+srcri^y^p,,)] 1 ' 

The investment cost is determined at the beginning of each 
stage and is given by the cost of altering network branches 
(by changing cross-sections of existing feeder sections or by 
installing new feeder sections), network nodes (by increasing 
the capacity of existing substations or by installing new ones) 
and charging stations. The cost of operation is considered at 
the beginning of each stage period and corresponds to the 
annual operational and maintenance cost of network branches 
that are in use, to charging stations, and to power that is 
curtailed. 

The total investment cost c' t nv (x, ,P t ,P t _ { ) at stage t is the 
linear function of decision variable x, , charging power P t 
and the charging power P t y (at stage t-l). It is calculated by 

inv / n n \ X 1 / s~iSQ SO , X 1 /~>SL SL \ , 

c, (x,,P t ,P l _ l )= 2^(C, x lt + 2^C, x,,) + 

I X X Z« + (2) 

I \ 1 (f-ip CO ^expand expand , fip (p _ p \ \ 

T Z-l fnetl P,' P P,' ^^ J vab\ r p,l r p,t-\>> 

pey/ c 

The total investment c" w (x t ,P t , P t _ { ) at stage t includes 
the substation capacity expansion fixed cost ( Cf°xf° ), and 

variable cost ( Cf L xf L , ) corresponding to choose candidate 
capacity of type SL, the capacity expansion cost of line j 
( C^x^ ) corresponding to choose candidate line type RJ, 

investment cost of building line k ( C^xff ) using line type 
AK, the fixed cost of building or expanding charging station 
at node p ( q ixed x c p 0 t , C e p x " a " d ' X ; p t and , respectively), and the 
variable cost of charging station investment 
( Cyab (P p t ~ P p i-\ ) )> which is proportional to the added 
charging capacity P pf -P pf _ x ■ 

The total operational cost c° per (r t , y t , P t , ) at stage t is 
the linear function of y t , charging station power P t and 
curtailed load power r t . It is calculated by the following 
function and the meaning of each part at the right side of 
equation is obvious. 

c°r{r t ,y t ,p t) = £ (°j°yj> X °fy^ 

j^W K J e ¥j 
key/ A ley/f P^¥° mey/ D 



(7) 



(8) 



C. Constraints 

The following constraints are considered: 

(1) Kirchhoff s laws, including KCL and Kirchhoff's voltages 
law (KVL). To obtain linear KVL constraints, the big number 
M is introduced to form the disjunctive version that is widely 
used in transmission expansion planning. The derivation of 
KVL can be found in [6]. 

s F f F + s R f R + s A f K +rt+gt = dt+ p t We 1; . . , J (4) 

\Zffr+[S F J V t \<M(l-yf t ) {ie ¥ F ,tel,-,T} (5) 

L J rowi ' 1 J 

| Z?f™ +[S* L V t |< M(l - y% ) {je ¥ R ,tel,-T} (6) 
\Zff^ + [s R l m V,\<M{\-yf tl ) 
{je ¥ R ,Je ¥ f ',fel,-7\} 
|Zf7*f K|<M(l-j<f) 
{ke V A ,Ke W t,te\,---T} 

(2) Operational limits 

• The constraints of maximum branch flow 

|^|<^/^,{Vz e ^,^l,..T.} (9) 

\ff, £yS/Zj+ Z >-".C,{v/ e rVe (io) 

I/m I* Z ygf^We ¥\te 1,-T] (11) 

• Constraints of substation capacity 

0<gf <gl_,+ Z (Z^kl,/ {V/e^,?el,-r.}(12) 

• Nodal voltage limits 

V^<V U <V^ WlG ¥ s ,tsl,-T.} (13) 

V m ^ m <V mJ <V imm {Mme ¥ D ,te\,-T.} (14) 

V^V ptt <V^ p {Vpe V c ,tel,-T.} (15) 

(3) Load shedding limits 

0<r ml <d mJ {Vms ¥ D ,tel,--T.} (16) 

(4) Charging station capacity limits, which set the limits of 
charging station capacity and total charging requirement. 

Z^^wrf, tel,-T. (17) 

^W^,.,^™,,, {V P e ¥ c ,tel,---T] (18) 

(5) Logic constraints, including all the logic constraints 
associated with decision variables. 

For each line that is candidate for replacement, only one 
type of conductor can be selected at all stages: 

ZZ^^W^^} (i9) 

<=1 Jsl/ff 

For each candidate line for addition, one type of 
conductor can be chosen at all stages: 



ZZ<f *l{Vker*} 



(20) 



1=1 Keii 



At any stage, only one capacity can be selected for 
substation expansion 

Z<f<l {Vle ¥ s ,tel,-T] (21) 

Lei//? 

No investment cost will occur if no candidate capacity is 
chosen: 

Z< <x?°{Vle ¥ s ,tel,-T) (22) 

Lswf 

For a line that can be replaced, only the selected 
conductor type can be used and operational costs will arise. If 
no candidate conductor is selected, the existing conductor 
will be used and result in operational costs. 

y% < x% jv/e ¥ R ,J e ¥ R ,te T} (23) 



^<i-X^{v/e^,^i,-r.} 



(24) 



At any stage t, only the selected conductor type can be 
used and operational costs will arise for the candidate new 
lines. 

yi K , < <f {\/ke ¥ A ,Ke ¥ i,te\,- -T] (25) 

If no charge capacity is added to the node p without 
charging facility at the initial stage, the charge capacity at this 
node should be equal to zero. 

I P PJ I^Z X Z ,{Yp e j/ 7 ,* e 1, ■ ■ T) (26) 

t=o 

For any charging node, the charging station can only be 
built once. If charging station is already there at the initial 
stage, the corresponding decision variable should be equal to 
0 at all the stages. 

± X c p °,<l{Vpe ¥ c } (27) 

(=0 

For a charging node p, building and expansion of the 
charging station cannot occur simultaneously. 

x expand + ^CO < jjy^ £ ^ £ j . . j, j (2g) 

At stage t, expansion is allowed when a charging station 
has been built: 

KT d * Z X Z {^pe ¥ c ,tel,- T) (29), 

where the x c p \ means if a charging station exists at initial 

stage when r = 0 . 

The capacity of the charging station is non-decreasing. 

0 < P Pit -P Pit _, < M(x c p 0 1+X ;'r d ){ype ¥ c ,tel,--T] (30) 

(5) Radial structure of distribution network. Improvements 
are made based on the work presented in [6] and the 
following constraints are derived. 

Z>- <U 'I Vtel-T (31) 



</ <X.v. Vtel,-T 

k<Ef A Kef( isy/ F jelff" Jey/f 

Z ft, Vf G l,...7\ 



(32) 



(33) 



The whole proposed formulation is a mixed-integer linear 
programming problem. It is solved by optimization software 
CPLEX®. 

IV. Simulation Results 

A. Test System 

The distribution network diagram and data can be found in 
[7]. The nodes 7, 10 and 13 are set as candidate places for 
building EV charging stations. All the data for the other 
nodes is the same as that given in [7], which is shown in 
Table I. The rated voltage of this system is 13.8kV and the 
allowable voltage range of each node is within [95%, 105%] 
of the rate voltage magnitude. 

The planning horizon is four years, divided into three 
stages, the first two being one-year duration and the third 

being for two years. S[" v = S° per =1, 8™ = 8 oper =0.9091, 
<5f v =0.8264, 8l per =1.5778. The total charging capacity 

requirements are 1.2 MVA, 6.0 MVA and 8.4 MVA for each 
stage, respectively. Parameters of candidate charging stations 
are listed in Table II. 



TABLE I 
Load Data of the Test System 



Node number 1 2 3 4 5 6 7 



9 10 11 12 13 14 15 16 



. , Stage 1 1.2 0 
Load 0i " , . „ , , 



/MVA 



0 1.2 1.2 1.2 / 1.2 1.2 / 1.2 1.2 / 0 0 1.2 

Stage 2 1.2 1.2 1.2 1.2 1.2 1.2 / 1.2 1.2 / 2.4 1.2 / 1.2 1.2 1.2 

Stage 3 1.2 1.2 1.2 1.2 1.2 1.2 / 1.2 2.4 / 3.6 1.2 / 2.4 2.4 1.2 

TABLE II 

Parameters of Candidate Charging Stations 



Investment 



Node 

number Fixed cost Fixed cost Variable 
(new)/ 1 0 3 $ (expansion)/ 1 0 3 $ cost/1 0 3 $/A 



_ Operational 



Capacity/MVA 



Cost/10 3 $/A Min Max 



7 


500 


100 


1 


1 


0.5 


2 


10 


500 


100 


1 


1 


1 


5 


13 


500 


100 


1 


1 


1 


5 



B. Results 

The planning results are illustrated in Table III and Fig. 1 . 
Totally, eight lines are built and three existing lines are 
replaced. The resulting network topology at each stage is 
strictly radial. 

The demands and charging requirements are all satisfied. 
Only two charging stations are built to reduce the investment 
and operational costs. One charging station is built at stage 1 
and expanded at stage 3. The other charging station is built at 
stage 2 and keeps unchanged at stage 3. 



TABLE III 

Summary of the Proposed Planning Results 



Stage 


Line 


Line 


Costs 
for 


Charging capacity (MVA) 


Charging 


Total 


built replaced network Node 

(10 3 $) 7 Nodel ° 


Node 
13 


costs 
(10 3 S) 


cost 
(10 3 S) 


1 


4 


2 


993 


0 1.5/new 


0 


717 




2 


4 


1 


646 


0 1.5 


4.5/new 


1146 


4713 


3 


0 


0 


24 


0 3.9/expansion 


4.5 


1187 





Q 0 0 0 



& 
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(a) Stage 1 



0 
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4.5MVA 
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(b) Stage 2 

— © © © 



0 



© 



(c) Stage 3 

Fig. 1 . Single line diagram and charging station capacity of the solution at 
each stage 

In order to show the merits of the multistage expansion 
method, the step by step planning method is simulated. The 



planning result obtained at stage t-l is set as the initial state 
of stage t. The corresponding planning results are given in 
Table IV and Fig. 2. Comparing the results obtained by the 
proposed method, the following differences should be noted 
although the final network topology is the same: 

(1) The step by step method builds and replaces more lines 
than the proposed method (11 and six, respectively) because 
the planning results in each stage are not co-optimized. The 
line between node 5 and 10 added in stage 1 is discarded in 
the following stages. Two lines added in stage 2 (line 10-11 
and 15-16) are not used in the final stage. 

(2) The capacities of the resulting charging stations are 
different. The charging station at node 10 is expanded twice. 

(3) The total cost of the planning result obtained by the 
proposed method is about 18% cheaper than that of the step 
by step planning method. 

TABLE IV 

Results Using Step by Step Planning Method 



q- 



Stage 


Line 


Line 


Costs for 
network 
(10 3 $) 


Charging capacity 
(MVA) 


Charging 
station costs 
(10 3 $) 


Total 
cost 
(10 3 $) 


built replaced 


N ° de Node 10 Node 13 


1 


5 


2 


609 


0 121 0 
new 


674 




2 


4 


1 


836 


1 92/ 

0 . 4.08/new 
expansion 


1336 


5555 


3 


2 


3 


471 


4 32 

0 4.08 
expansion 


1629 
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(b) Stage 2 
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(c) Stage 3 

Fig. 2. Single line diagram and charging station capacity of the solution at 
the each stage using step by step planning method 

In the above simulations, the total costs of 
building/expanding charging stations are higher than that of 
network expansion costs. Although the test data is not from 
practical projects, the costs of building charging facilities and 
the role of coordinated planning between distribution 
networks and charging facilities should not be neglected. 

V. Conclusions 

A mathematical formulation is proposed in this paper to 
solve distribution network expansion planning and optimal 
siting and sizing of EV charging station at the same time. A 
mixed-integer linear programming problem is built with 
comprehensive constraints considered, which includes the 
constraints of voltage magnitude limits and radial network 
topology requirement. 

Simulations are carried out on the 18-node test system. Test 
results show that the proposed method can reduce the grid- 
connection cost of satisfying EV charging requirements. 
Comparisons show the superiority of the proposed multistage 
expansion planning method. More work is needed to put the 
method into practical use. 
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